Experimental demonstration of bosonic commutation relation via superpositions of 

quantum operations on thermal light fields 
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We present the experimental realization of a scheme, based on single-photon interference, for 
implementing superpositions of distinct quantum operations. Its application to a thermal light field 
(a well-categorized classical entity) illustrates quantum superposition from a new standpoint and 
provides a direct and quantitative verification of the bosonic commutation relation between creation 
and annihilation operators. By shifting the focus towards operator superpositions, this result opens 
interesting alternative perspectives for manipulating quantum states. 
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The superposition principle is one of the pillars upon 
which the entire structure of quantum mechanics is 
built [H. A quantum system in a pure state can always 
be described as a superposition of linearly independent 
states; thus once one has a quantum system represented 
by a pure state, the superposition is naturally there. An 
inter-body superposition state, the so-called entangled 
state, is somewhat trickier to generate than a single-body 
superposition state. However, it has been demonstrated 
that entanglement can be achieved by various methods, 
including a series of unitary operations [3, IJlor by 
post-selection of events after unitary operations [5| . On 
the other hand, the discussion about superpositions of 
classical mixed states is not as clear as for a pure state Q . 

Quantum operators, besides quantum states, play a 
crucial role in describing physical operations including 
unitary transformations and measurements in quantum 
theory. If one can implement a superposition of opera- 
tors, one can also construct state superpositions by ap- 
plying the superposed operators to a given state, unless 
it is a simultaneous eigenstate of the component opera- 
tions. In fact, also the Schrodinger's cat paradox [7| can 
be understood as the quantum-mechanical impact of the 
superposition of macroscopically distinct operations (to 
kill or not to kill) on a classical object (the cat). 

Several groups have recently succeeded in applying 
simple quantum operators to different quantum states. 
For example, in the optical domain, basic operations, 
like single-photon addition and subtraction have been 
demonstrated to produce highly nonclassical 
and non-Gaussian states 12] even when applied to clas- 



detection of a single photon in an ancillary (herald) light 
mode. Sequences of photon additions and subtractions 
have also been implemented to show that the two se- 



quences aa 



^ and 



where and a are the bosonic 



creation and annihilation operators, give different results 
when applied to the same input light state 15|. This is 



an important corner stone for the proof of the bosonic 
commutation relation 



[a, a^] = aa^ — a^a = 1 



(1) 



sical states of light 13, 14l|- Both photon addition and 
subtraction are performed in a conditional way upon the 



which is at the heart of many important consequences 
of quantum mechanics. However, the complete demon- 
stration of the commutation relation was out of reach 
because of the lack of an important element in the quan- 
tum manipulation toolbox: the possibility of superposing 
different operators A and B to form the general opera- 
tor superposition aA + (3B, where a and /3 are complex 
amplitudes. 

Since the superposition principle relies on the indistin- 
guishability among different alternatives, the experimen- 
tal implementation of quantum operators heralded by a 
single-photon detection offers a very convenient way to 
achieve this goal. If the herald field modes of different 
operators are properly mixed by means of a beam split- 
ter, the information about the origin of a click in the 
herald photodetector is erased and a coherent superposi- 
tion of the different operators can be conditionally imple- 
mented. Somewhat similar schemes have been recently 
proposed and experimentally implemented for increasing 
the entanglement of bipartite Gaussian quantum states 
by inconclusive photon subtraction 16, 1^ [l^, toward 
the implementation of a quantum repeater for l ong -haul 
quantum communication in ionic systems [20f|, and 
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for the remote delocalization of a single photon over dis- 
tinct temporal modes 



21| 



In this Letter we present the experimental realization 
of a general scheme, based on single-photon interference, 
for superposing distinct quantum operations. As recently 
proposed in Ref. demonstrating the bosonic com- 

mutation relation thus reduces to realizing the balanced 
superposition aa^ — d^d, and showing that it corresponds 
to the identity operator 1. While most of the mathe- 
matical structure of quantum mechanics is based on the 
commutation relation, this is the first time it is directly 
probed in an experiment. 

The primary laser source is a mode-locked Ti:Sa laser 
emitting 1.5 ps pulses at a repetition rate of 82 MHz. A 
rotating ground glass disk (RD) is inserted in the path 
of the laser beam and a bare single-mode fibre (SMF) 
is used to collect a portion of the scattered light to pro- 
vide the pulsed thermal light states of mean photon 
numbers around unity which have been used as the input 
states in the experiment. We use a convenient modular 
scheme (see Figure 1), where one single-photon addition 
stage 'lO|| (at, based on conditional single-photon para- 
metric amplification in a type-I BBO -(3 barium borate- 
nonlinear crystal) is placed between two single-photon 
subtraction stages (a, based on the conditional re- 
moval of a single photon by adjustable low-reflectivity 
beam splitters, obtained with combinations of half-wave 
plates (HWP) and polarizing beam-splitters (PBS)). 

The successful implementation of the desired superpo- 
sition of operator sequences is determined by the coinci- 
dence (C) between a click from the single-photon detec- 
tor heralding photon addition (Da) and one from a sec- 
ond photodetector (Ds), placed after a -3db fibre coupler 
(FC) combining the herald signals from the two subtrac- 
tion modules. By a click from detector Ds we know that 
a single photon has been subtracted, but we cannot iden- 
tify if it was before or after the photon addition. In these 
conditions, a coincidence event heralds the application of 
the general operator superposition with an 

adjustable phase 0, to any input light field. By varying 
the phase (j) with a piezo-actuated mirror, any arbitrary 
balanced superposition of the two operator sequences can 
be obtained. In particular, by setting (j) — or (j) — n, 
one can directly implement the commutator or the anti- 
commutator of the creation and annihilation operators, 
respectively. Note that the present scheme, differently 
from the original theoretical proposal [2^ . where only 
these two possibilities were allowed by using two detec- 
tors at both exits of the beam splitter, allows for greater 
flexibility (by generating operator superpositions with ar- 
bitrary relative phases) and experimental simplicity. 

Without the click from the addition module (a^), the 
scheme reduces to a Mach-Zehnder-type interferometer 
which can be used to verify the indistinguishability of 
the two subtraction events by evaluating the visibility of 
the single-photon interference in the counts at detector 
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FIG. 1: Experimental setup. Pump pulses for paramet- 
ric down-conversion are obtained by frequency-doubling the 
laser output in a LBO -lithium triborate- crystal. The 
conditionally-prepared signal state is mixed with a strong ref- 
erence coherent field (LO, obtained from a portion of the main 
laser output) on a 50-50 beam splitter whose outputs are de- 
tected by two photodiodes (Hamamatsu PIN S3883). The 
balanced homodyne detection (BHD) signal is acquired and 
stored by a digital oscilloscope (Tektronix TDS7104) on a 
pulse-to-pulse basis triggered by a coincidence (C) between 
clicks from the Da and Ds single-photon detectors (Perkin 
Elmer model SPCM AQR-14). 



Ds. Visibilities of about 97% are obtained by carefully 
balancing the reflectivity (« 3%) of the involved beam 
splitters, by fine polarization control (PC), and by a pre- 
cise adjustment of the delays between the corresponding 
herald modes. If delays are not compensated, only a sta- 
tistical mixture of the two operator sequences with equal 
weights is obtained. The Ds count rate is also used to 
monitor the superposition phase (j) a-nd lock it to any de- 
sired value. The effects of experimental deviations from 
the ideal realization of this scheme (such as the finite re- 
flectivity of the subtraction beam splitters, the possible 
multiple photon-pair production in the parametric pro- 
cess, or the fact that real photodetectors are not able 
to discern the number of photons but only there being 

not to 



mi 



photons or not) have already been shown 
significantly affect the results of the experiment for the 
present range of parameters. 

The state resulting from the chosen operator superpo- 
sition is analyzed by means of a high-frequency, time- 
domain, balanced homodyne detector [2^ yielding the 
distributions of measured field quadratures. Since both 
the initial thermal states and those resulting from the 



above manipulations possess no intrinsic phase, the phase 
of the local oscillator (LO, the reference coherent field for 
homodyne detection) is not actively scanned, and phase- 
independent marginal distributions are obtained. How- 
ever, the final states still clearly depend on the phase 
(j) of the superposition. High experimental efficiency is 
obtained by minimizing all spurious losses and making 
sure that all the single-photon operations are performed 
in exactly the same spatiotemporal mode as the one se- 
lected by the LO. This requires narrow spatial and spec- 
tral filtering (F) in the herald mode of the parametric 
down-conversion crystal, and an accurate matching of 
the fibre-coupled fields reflected from the two subtracting 
beam splitters to the LO spatial mode. 

Figure 2a) shows a sequence of histograms of raw ho- 
modyne data acquired while scanning the phase of the 
superposition. The remote manipulation of the state 
by the implementation of different superpositions of cre- 
ation and annihilation sequences is clearly observed. The 
quadrature distribution of the final state undergoes a 
very rapid initial evolution from a bell-shaped curve at 
(/i = towards a volcano-shaped one around (j) — n, where 
the phase dependence is much slower. The phase change 
would not have resulted in such different output states if 
the operations had been statistical mixtures. 

The special cases of </> = and </) = tt are illustrated 
in more detail in Figure 3, where the Wigner functions 
of the original thermal state and those resulting from 
the experimental realization of the commutator and anti- 
commutator between a and are presented. The fi- 
delity F = \Tr^/y/pi^poutVpin\'^ between the original 
thermal state and the final one (represented by their re- 
constructed density operators pin and Pout, respectively) 
is about F = 0.992 for the commutator case (0 = 0). 
This demonstrates that the implemented operator su- 
perposition is essentially equivalent to the identity op- 
erator. Wigner functions have been obtained from the 
fO diagonal density matrix elements (13 for the anti- 
commutator case) reconstructed by means of a maximum 
likelihood algorithm without any correction for 

the finite detection efficiency. If homodyne detection effi- 
ciency {rjd — 0.7) is corrected for, the Wigner function for 
the state resulting from the anti-commutation operator 
{(j) = it) clearly attains negative values. The existence of 
negative regions in the reconstructed Wigner function is 
a direct signature of the fact that the state impinging on 
the homodyne detector is highly nonclassical. 

Actually, it is interesting to note a fact that was 
not realized in the theoretical proposal [i^]: because 
of the normalization of quantum states, the above re- 
sults just demonstrate the commutation relation up to 
a multiplicative constant K, i.e., one might still have 
[a, a^] = Kl. However, in this case, the anti-commutator 
setup implements the 2a^a -I- Kl operator, which pro- 
duces an output state strongly depending on the exact 
value of the constant K. Figure 2b) reports the mea- 
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FIG. 2: a) Sequence of histograms corresponding to about 
250,000 raw quadrature data for the final state as a function 
of the normalized quadrature x and of the superposition phase 
<j!>. Data have been obtained in a 16-hour measurement, and 
have been binned in 9 phase intervals between and tt. A 
thermal field with a mean photon number n — 0.9(1) is used 
as the initial state, b) Histogram of raw quadrature data 
(solid dots) for the anti-commutator setup at (j) — tt. Also 
shown are theoretical curves calculated for the actual exper- 
imental parameters (total efliciency rj = 0.61, n = 0.9) and 
different values of the commutator {K — 0: dashed orange; 
K = 2: dotted green; K = dash-dotted blue). The solid 
red curve is the result of the best fit to the experimental data. 



sured homodyne quadrature distribution for the same 
initial thermal state after the application of the anti- 
commutation operator {(j) — tt). Also reported are the 
theoretical distributions calculated for the same exper- 
imental parameters but with a few different values of 
the constant K . Experimental data are consistent with 
K = 1^ whereas different integer values are in evident 
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FIG. 3: Experimentally reconstructed Wigner functions of the 
original thermal state and of those resulting from the applica- 
tion of the commutator and anti-commutator superpositions. 
The state is not changed for (j) — 0. About 10* (10^) quadra- 
ture data points have been acquired in the (anti-)commutator 
case. 

disagreement. A best fit of the experimental homodyne 
data gives K = 1.02(3), thus quantitatively demonstrat- 
ing the bosonic commutation relation. 

Although the present case only required a coherent 
superposition of two (sequences of) quantum operators 
with the same weight, the proposed scheme is much more 
general and allows one, in principle, to implement co- 
herent superpositions of an arbitrary number of opera- 
tors with arbitrary relative amplitudes and phases by a 
network of beam splitters with adjustable reflectivities. 
The single-photon interference as a way to produce gen- 
eral operator superpositions can be straightforwardly ex- 
tended to the superposition of a and by letting the 
creation and annihilation herald photons interfere at a 
beam splitter. As a representative example, "position" 
and "momentum" operators in the phase space can be 
implemented in this way. Our approach can even be gen- 
eralized to realize various superpositions of higher-order 
operators in terms of a and . Any quantum state can 
be written as C„a^" |0), where C„ are complex ampli- 
tudes, and arbitrary states can thus be generated by ap- 
plying the appropriate superposition of photon creation 
operators. In this paper we have experimentally demon- 
strated a basic building unit for such general operator 
superpositions on a traveling light field for the first time. 
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